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ABSTRACT 

AGN outflows are the heat given up when gas in a galaxy evolves towards ther- 
modynamic equilibrium. Indeed, while AGN feedback regulates the growth of massive 
galaxies, its origins can be understood as the spontaneous thermodynamic process 
which ensures that the (Gibbs) free energy of the system always decreases, enabling 
the galaxy to reach a more energetically favourable state. In particular, it is shown 
that feedback heating processes will be favoured whenever the hot atmosphere of a 
galaxy would effectively gain energy as a result of cooling. For example, as the hot 
atmosphere of a galaxy cools and contracts, the work done by gravity will be ther- 
malised, with a fraction of the gas also being captured by stars and the supermassive 
black hole at the centre of the galaxy. If this gain of energy exceeds the loss of energy 
that occurs when cooling gas drops out of the atmosphere, the Gibbs free energy of the 
system would increase overall. Since this is energetically unfavourable, feedback heat- 
ing is initiated which acts to reduce the net cooling rate of the atmosphere, thereby 
preventing any build-up of energy. The Gibbs free energy can also decrease in the 
absence of feedback heating, but only if the loss of energy due to mass dropping out of 
the atmosphere exceeds the gains of energy described above. Therefore, to ensure that 
the Gibbs free energy always decreases, a galaxy will necessarily flip between these 
two states, experiencing episodes of heating and cooling. Due to the close long-term 
balance between heating and cooling, the gas in a galaxy will evolve quasi-statically 
towards thermodynamic equilibrium, which has the observable appearance of galaxy 
growth being regulated by AGN feedback. The same mechanism also provides an ex- 
planation for why strong AGN feedback occurs more frequently in cool-core galaxy 
clusters than in non cool-core clusters. 
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1 INTRODUCTION 

The cooling time of X-ray emitting gas near the centres of 
many massive galaxies and galaxy clusters is much shorter 
than the Hubble time. In the absence of heat sources, sig- 
nificant quantities of the gas would cool and form stars. 
However, X-ray spectroscopy has shown that the rate at 
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that the gas is somehow being reheated. 

A wealth of both observational and theoretical evidence 
suggest that energy input by a central Active Galactic Nu- 
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cleus (AGN) is responsible for reheating the gas. In particu- 
lar, powerful outfiows produced by AGN interact strongly 
with their environment, inflating lobes of radio emission 
in the hot gas that permeates massive g alaxies and clus- 
ters of galaxies ( e. g. Birzan et 200 j : Best et all 
Dunn et al.l I2OO5I: [Raffertv et al.l ]2006l : [ Best et al 



2005; 



20071 : 



Shabala et al.ll2008l ). In addition, recent observational stud- 
ies of Brightest Cluster Galaxies (BCGs) suggest that this 
radio AGN activity is related to the thermal state of its en- 
vironment. Systems with short radiative cooling times, or a 
low central entropy, are more likely to exhibit active star for- 
mati on, opt i cal line-emission and o utfiow-producing AGN 
(e.g. Bur nd Il99d: ICrawfor d et~al] Il999l: [Cavagnolo et all 
2008; Mit tal et al.| [2009: Ra ffertv et al.ll2008h " 



Imple menting and extendin g the ideas of early 
work (e.g. iBinnev fc TaboJ Il995l : iTucker fc DavidI 119971 : 
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ICiotti fc Ostrikeij|200ll : ISilk fc R ees"l998') , theoretical stud- 
ies have drawn attention to the potential, wid e-ranging im- 
pact o f AGN feedback (e.g . [Pabian 1999; Briig gen fc Kaisei 
20021 : lOmma et al.l 12004: iKina i2005i: iSiiacki fc Springe] 



200e : iDi Matteo et al.l l2008l : iBooth fc Schavel l2009l : ISoker 
2OO9I : iMcCarthv et al.ll201Gl : iGaspari et al.ll201ll ). Notably, 



semi-analytic models of galaxy formation have demon- 
strated that, in principle, AGN heating can both re- 
heat cooling flows and explain the exponential cutoff at 
the bright end of the galaxy luminosity function (e.i 
iBenson et alll2003l: ICro ton ct al.l l2006l : iBower et aLll2006 
see also l|Short fc Thomas, , 20091 ). In addition, AGN heat- 
ing has been shown to be crucial in shaping the X-ray 
luminosity-tem perature rel ation of mass i ve ga l axies (e.g. 
iPuchwein eTa l. 2008; Bow er et al] l2008l : iPopel |2009| ). As 
such, AGN feedback has become a ubiquitous feature of 
galaxy evolution models, underpinning much of our under- 
standing of why massive galaxies are they way they are. 

However, the accuracy of galaxy evolution models is 
limited by their spatial resolution and our incomplete un- 
derstanding of important physical processes - particularly 
those that occur on sub-grid spatial scales, such as the the 
microphysics of the gas, and black hole accretion. As a result, 
it remains difficult to implement feedback in a fully self- 
consistent manner, hindering our understanding of galaxy 
evolution. For this reason, it is vitally important to better 
understand the fundamental processes that drive and trigger 
feedback. 

As an example of an AGN-triggering process that is not 
curre ntly captured in c osmological models of galaxy forma- 
tion, iPope et all l|201ll ) demonstrated that gas in the cen- 
tre of a massive galaxy can become gravitationally unstable 
on sub-kiloparsec scales, thereb y providing the -AGN w ith 
an abundant source of fuel (see ISokeil 120061 : iKindlioO^ . for 
different examples) . The chara c terist ic duty cycles AGN ac- 
tivity predicted by Pope et al. (I2OII ) are in agreement with 
obser vations (e.g. lBest et al.ll2005l : lBestll2007l : ishabala et all 
|2008| ). suggesting the model does describe some important 
features of the AGN feedback phenomenon. As a result, this 
knowledge can be used to develop an improved prescription 
for AGN activity in a numerical model that does not explic- 
itly account for self-gravity on sub-kiloparsec scales. 

If correct. Pope e t al.. (^2011,) goes some way to explain- 
ing an important physical process that triggers AGN feed- 
back in a hot atmosphere, but it does not explain why the 
system should evolve in such a way: why does AGN feedback 
occur at all? 

To re-phrase the problem: assuming that a galaxy 
can be considered to be a thermodynamic system, it 
should naturally e volve to minimis e its free ene rgy (e.g. 
iLvnden-Bell fc Wo od 19680 (see also lNusserll2009l . for a dis- 
cussion of galaxy cluster thermodynamics without feedback 
heating). With this in mind, it is perhaps surprising that 
AGN feedback occurs at all since, by re-heating the galaxy's 
atmosphere, it appears to be obstructing the galaxy's ten- 



^ More specifically, a galaxy is likely to be an open thermody- 
namic system since it can exchange mass and energy with its 
environment. In contrast, a closed thermodynamic system is one 
that only exchanges energy with its environment, while an iso- 
lated system shares neither mass nor energy with its environment. 



dency to evolve towards a minimum free energy configura- 
tion. 

At this point, it is necessary to distinguish between 
AGN heating which is driven by external inffuences such 
as galaxy mergers, and AGN heating that is driven inter- 
nally by various properties of the host galaxy such as the 
X-ray luminosity, temperature and mass of its gaseous at- 
mosphere. The latter constitutes AGN feedback in the truest 
sense of the definition and will be the focus of this paper. 
The distinction is relevant because merger events push the 
system further away from thermodynamic equilibrium. In 
contrast, internally-driven processes occur spontaneously as 
the system evolves towards thermodynamic equilibrium by 
reducing its free energy (e.g. lFermilll956l : lLandau fc Lifschitj 
|l980). However, it is important to mention that merger- 
driven AGN heating may also play a significant role in 
galaxy evolution, with recent work suggesting that the ef- 
fect may be more conspicuous for galaxies in group and field 
environments than in clusters (see Kavirai et al.ll201ll . for a 
recent example). 

As indicated above, the progression towards thermody- 
namic equilibrium provides an informative framework with 
which to describe the characteristic features of gas evolution 
in a galaxy. Indeed, it is important to note that while the 
atmosphere of a galaxy may be close to hydrostatic equilib- 
rium it is far from thermodynamic equilibrium either with 
its environment or internally. This is because the atmosphere 
is much hotter than the Cosmic Microwave Background, 
and it contains stars and a super massive black ho le which 
have a range of temperatures (e.g. lHawkin3ll976l ). There- 
fore, according to the second law of thermodynamics, there 
will be large fluxes of heat from the hot atmosphere to its 
environment, and to the stars and the supermassive black 
hole. In particular, the very low temperature of a super- 
massive black hole means that any accretion of mass and 
energy from its environraent results in a very large entropy 
increase (e.g. lBekensteinlll973l : lBardeen et alll 19731 : iHawkind 
Il976i ) and is, thus, highly probable. However, the power 
output and subsequent heating that arises when the black 
hole grows prevent the accretion rate from rising uncontrol- 
lably. The competing influences of energy and entropy are, 
most generally, quantifl ed by the free energ y of the system 
(e.g. ,Ferml,1956i : iLandau fc Lifsch"it3ll980l ). Consequently, 
the free energy can be considered to be a thermodynamic 
potential. To date, this approach has proved to be useful 
in describing stable states o f gravitating g as spheres (e.g. 
Lvndcn-B eU fc Woo d 1968; Stabler 1983; Tomisak a et al.l 
.1988,: McKee fc Hollimaa,1999 : Bludman fc Kennedvj201ll ). 
In this article, we flnd the Gibbs free energy to be the most 
suitable form of thermodynamic potential with which to de- 
scribe and explain the driving force behind the feedback 
processes thought to be key in governing the evolution of 
the galaxy. 

Formally, a thermodynamic system will achieve a state 
of thermodynamic equilibrium when its free energy reaches 
a minimum (e.g. iFermi 1955 ). As a result, for systems that 
have not yet reached equilibrium, processes that reduce the 
free energy will be favoured and will occur s pontaneously 
(e.g. iFermil [l956l : iLvnden-Befl fc WoodI [l968l ) . In this for- 
malism, AGN and stellar feedback are the spontaneous pro- 
cesses which ensure that the free energy of the system always 
decreases, so that the system tends towards thermodynamic 
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equilibrium. Therefore, it is shown that AGN feedback is 
energetically favourable. This is primarily because feedback 
heating prevents the effective build up of energy that occurs 
as hot gas from the galajcy's atmosphere cools and flows to- 
wards the centre of the galaxy. As this happens, the stored 
gravitational potential energy of the atmosphere is converted 
into heat, while significant quantities of the cooling gas will 
be gravitationally captured by stars and the supermassive 
black hole at the centre of the galaxy. Gas that is gravita- 
tionally captured by a compact object will be bound to it 
with an energy that is a non-negligible fraction of its rest 
mass energy. Since this is much greater than its binding en- 
ergy to the host galaxy, the gas can effectively gain energy as 
it cool^. Finally, as the gas falls onto the compact object, 
the binding energy will be radiated as heat, reducing the 
net cooling rate of the atmosphere. Consequently, feedback 
prevents the build up of energy that can occur when gas 
cools in the presence of compact objects, and ensures that 
the galaxy evolves slowly towards a minimum free energy 
configuration. 

The article is organised as follows: in section 2 we derive 
a general thermodynamic galaxy model, in section 3 we con- 
sider the broader implications of the model, and summarise 
the findings in section 4. 



2 MODEL 

While this model mainly focuses on explaining the origin of 
feedback heating in massive elliptical galaxies, the principles 
also hold for stellar feedback, and AGN feedback in galaxy 
clusters. 

To ensure the study is applicable to the widest range 
of environments we start with a thermodynamic description 
of a hot galaxy atmosphere (with an embedded cool phase) 
that flows within the galaxy's darkmatter gravitational po- 
tential. The gravitational field exerted by the darkmatter 
is comparatively smooth and extends far beyond the visible 
extent of the galaxy, ft is this field that governs the large- 
scale properties of the hot atmosphere. Differences in the 
global heating and cooling rates drive the transfer of mass 
between the hot and cool phases. Gas belonging to the cool 
phase is assumed to be strongly influenced by the presence of 
stars and the black hole since it can be more easily accreted 
than hotter gas, due to its lower thermal velocity. For this 
reason, the approach is extended to consider a generalised 
galaxy model which describes the 3 main gravitational in- 
fluences on the gas: i) the darkmatter; ii) the stars; iii) the 
supermassive black hole. The system then favours heating 
and cooling process which cause the galaxy's gas to move 
between the gravitational components, re-distributing mass 
and energy, in such a way that it evolves towards a thermo- 
dynamic equilibrium state. 



^ This phenomenon is simil ar to th e nega tive heat capacity of 
galaxy clusters discussed by 'Nu sserl | |2009| ). However, it differs 
because feedback heating intervenes to prevent the heat capacity 
becoming negative. 



2.1 Total energy of a hot gaseous atmosphere 

The hot gaseous atmosphere that permeates a massive 
galaxy is thermally ionised, and cools through the emis- 
sion of X-rays. As the gas radiates energy it loses pres- 
sure support so that the weight of the overlaying gas 
causes an inflow. Nevertheless, the gas is generally as- 
sumed to be close to hydrostatic equilibrium because the 
coo ling time is much longer t han the free-fall time (e.g. 
see iMcNamara &: NulsenI |2007| . for a review) . Additional 
deviations from the hydrostatic limit are caused by bulk 
flows and turbulent eddies driven by m erger events and 
large-scale outflows fro m the galaxy (e.g. [R ickcr fc Sarazi 
, 2001; Diehl fc Statler 20o3: IScannapieco fc Briiggcn 2l 
Isternberg fc Sokeill2009l : iDong et al.ll2010l : iPellegrinilboill ) 
As such, the evolution of a hot atmosphere can be described 
in terms of its kinetic, thermal and gravitational potential 
energy. 

For the purposes of this study, it is necessary to define 
a region of interest within the galaxy atmosphere. The stan- 
dard choice is the volume enclosed by the cooling radius 
- the radius at which the gas cooling time equals (a frac- 
tion of) the Hubble time. The reason for this choice would 
be that the volume enclosed by the cooling radius provides 
an estimate of the gas mass that can have cooled within 
the lifetime of the Universe. The cooling radius, rcooi, often 
implicitly enters calculations which relate the bolometric X- 
ray luminosity of the hot gas, L, within rcooi, its average 
temperature, T, within the same volume, and the so-called 
classical mass flow r ate, ritciaa- Th e magnitude of this flow 
rate is given by (e.g. lFabiarj[l994l ) 



L 



7 fcbT 



(1) 



(7 - 1) (m) 

where 7 is the adiabatic index of the hot gas, fcb is Boltz- 
mann's constant and (m) is the mean mass per particle. 
According to the derivation of equation ([1} from the fluid 
energy equation, mdas is the mass flow rate due to the con- 
traction of the hot gas of temperature T as it radiates energy 
at a rate L. As such, rhdas is the rate at which mass flows 
through rcooi due to the contraction of the hot gas within 
^cooi ■ More generally, the mass flow rate across rcooi is deter- 
mined by the difference between the global heating rate, H, 
and the cooling rate, L, such that H — L oc mtot ■ In the limit 
that the change in gravitational potential energy is small by 
comparison, the mass flow rate of the hot atmosphere takes 
the following form (e.g. Sarazin 1986) 

1 



H ~L : 



-mtotCs 



(2) 



(7-1) 

where Cg = '^ki.T/ (m) is the sound speed of gas in the atmo- 
sphere. Thus, ii H < L, the atmosphere cools overall, and 
mass will flow inwards, li H > L, the gas is heated overall, 
and mass will flow outwards. It is important to note that 
H is the rate at which the injected energy is dissipated as 
heat; as such, it represents a fraction, rjc, of the total en- 
ergy injected to the system. For example, an AGN accreting 
material at a rate less than ~3% of the Eddington limit is 
likely to be radiatively inefficient so that the majority of 
the power output is in the form of kinetic outflows that are 
thought to couple strongly to the ambient gas, i.e. »7c ~ 1. 
Conversely, accretion above this critical rate is radiatively 
efficient, meaning that the power output is predominantly 
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in the form of photons. In the radiatively efficient mode, it 
has been argued that only 77c ~ 5% of the accretion power is 
avail a ble to heat the surrounding gas (e.g. ISiiacki fc Springel 
I2OO6I : iKind |2009| '1 . The coupling parameter rjc can also be 
used to account for the fraction of energy dissipated within 
the region of interest for outflows that extend beyond its 
outer boundary. 

While the cooling radius provides a reasonable choice 
for an outer boundary under many circumstances, the 
changing mass within rcooi and the unknown variation of the 
external pressure at rcooi make calculations of the work done 
on the region of interest much more complex than they need 
to be. For this reason, the calculations presented in this in- 
vestigation make use of a more appropriate outer boundary 
condition that: i) simplifies the calculation of work done on 
the region of interest; ii) ensures there is no mass flow across 
this boundary, as shown below, meaning that the system is 
formally closed. In particular, we choose an outer bound- 
ary which is deflned by a virtual surface at flxed pressure 
(Axt) that moves inwards and outwards depending if the gas 
within is cooled or heated overall. To prove that there should 
be no mass flow across an outer boundary of flxed pressure, 
for quasi-static changes, we present the following argument. 
For a region of gas bounded by an external pressure Pcxt, 
heated at a rate H and cooled at a rate L, the rate of change 
of volume is given hy H — L = PcxtV'7/ (7— !)■ Using the ideal 
gas law, the pressure is written Pcxt ~ pk\^T/{m), where p 
and T are the density and temperature at the boundary. 
Here, we assume that the temperature throughout the at- 
mosphere is constant, so that the temperature at the outer 
boundary is equal to the average temperature. Then, assum- 
ing that the region is spherical, the rate of volume change 
is 1/ = A-rrr^r, with r being the rate of change of radius at 
which the external pressure is Pcxt- Substituting for Pcxt and 
V gives 

rr r ,1 2 . 7 1.2 /ON 

(7-1) (m) (7-1) 

Thus, for an appropriate choice of Pcxt, the gas volume of 
interest will change radius at approximately the same ve- 
locity as the local mass flow given by equation 0, mean- 
ing that m ~ rhtot- Therefore, under ideal conditions, there 
should be no mass flow across the boundary deflned by con- 
stant pressure, but it is important to remember that, due 
to differences between the global heating and cooling rates, 
there is still mass flowing within the region, i.e. relative to 
a flxed radius. In particular, the mass flow rate relative to a 
flxed cooling radius has a magnitude of ~ mtot . For this rea- 
son, and to ensure the region is lar ge enough that feedback 
energy is dissipated with in it (e.g. iHeinz et al.l l2006l : iPopd 
I2OO9I : iFabian et al.|[201ol ). Pext is taken to be the pressure 
at the initial cooling radius. Finally, the region of interest 
is assumed to be spherical and centred on the supermassive 
black hole at the centre of the galaxy. Under these assump- 
tions, the total energy of the gaseous atmosphere within the 
region of interest, n, can be written as 

E{Q.) = K{Q.)+Q{Q.) + ^l!{Q.), (4) 

where K is the bulk kinetic energy, Q = U -\- PcxtV^ is the 
enthalpy of the gas volume of interest, with U being the 
internal energy and V being the volume of O, and ^ is the 
gravitational potential energy. 



As the gas within cools, the gas will lose gravitational 
potential energy while simultaneously having compression 
work done on it. Initially, both the gravitational and com- 
pression work will increase the kinetic energy of the gas. 
However, as long as the relaxation time of the gas is short 
compared to the work timescale, the kinetic energy will be 
thermalised within Q.. In this limit, the kinetic energy of the 
system will remain constant, the system will evolve quasi- 
statically, and the compression can be considered to be re- 
versible. The relationship between the kinetic, thermal and 
potential energies can be derived by differentiating equation 
Q with respect to time 

E^k + e + ii, (5) 

where we have omitted Q, for brevity. 

According to the argument above, the rate of change of 
kinetic energy is equal to the total rate at which work is done 
on the system, less the rate at which the kinetic energy is 
thermalised, D. The total rate at which work is done on the 
system is sum of the rate at which gravity does work (—3') 
and the rate of compressional work done on the system (VK), 
such that K — W — '^ — D. Then, in the limit that the work 
is performed reversibly within fi, we can write D = W — . 

By deflnition, the rate of change of enthalpy is given by 
the sum of the heat fluxes, Q, plus the kinetic energy dissipa- 
tion rate, less the rate of compression work; G) = Q + D — W. 
Thus, in the limit of reversibility, we have O — Q — ^P, which 
is the flrst law of thermodynamics for a system in which the 
change of gravitational potential energy is important. Fi- 
nally, the rate of change of total energy is E — Q + W. 

Under the assumption that the hot gaseous atmosphere 
evolves quasi-statically, it is appropriate to use further ther- 
modynamic arguments to describe the direction of processes 
driving this evolution. These arguments are outlined and ex- 
plored below. 

2.2 The Gibbs Free Energy of a galaxy 
atmosphere 

A system will be in stable thermod ynamic equi librium when 
its free energy is at a minimum re.g. lFermill95^ ). As a result, 
processes that reduce the free energy will be favoured and 
will occur spontaneously (e.g. iLvnden-Bell fc Wood 196^). 
In contrast, processes that increase the free energy are en- 
ergetically unfavourable and will be improbable. 

The appropriate form of the free energy for a system at 
constant pressure Pcxt and temperature T, is the Gibbs free 
energy - henceforth referred to as GFE - which is deflned 
as 

G^e-TS (6) 

where 9 is the enthalpy as deflned previously, T is the sys- 
tem temperature and S is the system entropy. 

Equation © should be viewed as a thermodynamic po- 
tential that quantifles the relative influences of enthalpy 
and entropy. Thus, for an isolated system (with constant 
U, Pcxt and V^), the GFE is minimised by maximising the 
entropy. In contrast, for a closed system that exchanges 
energy with its environment, the GFE will be minimised 
by de creasing the enthalpy and increasing the entropy (e.g. 
iLvnden-Bell fc Wood 1968). 

Following the deflnition of the GFE, a system will be 
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in equilibrium when G is a minimum; this occurs when 
the change in GFE is zero , i.e. AG = (e.g. Feriiii|ll956l : 
iLvnden-Bell fc Woodlll96^ 'l. Using equation ([5}, the change 
in the GFE is written 



AG = Ae - TAS - SAT. 



(7) 



To generate a more useful form of the differential GFE, we 
substitute for AO and TAS, and assume that the temper- 
ature of the atmosphere remains approximately constant so 
that SAT = 0. This assumption is justified because the cool- 
ing time of the gaseous atmosphere is generally much longer 
than the gravitational free-fall time. Therefore, whether the 
gas is being heated (expanding) or cooled (contracting), the 
gas will flow sufficiently slowly that it can be considered 
to be close to hydrostatic equilibrium, i.e. its sound speed 
must comparable to the velocity dispersion of the underly- 
ing gravitational potential. Since the gravitational potential 
is taken to be static, the gas temperature must also be ap- 
proximately constant. 

From its definition, the change of enthalpy is written 
AG = AU + PextAl/ + VAPcxt. Since the outer bound- 
ary of Q is defined to be movable so as to maintain con- 
stant pressure (APcxt = 0) the differential enthalpy becomes 
A0 — AU + Pcxt AV. Furthermore, the second law of ther- 
modynamics for the reversible transfer of heat, AQ, from a 
system at temperature, T, to its environment, is written 



AQ = TAS 



(8) 



where AS" is the entropy change of the system. Accordingly, 
equation ((TJ becomes 



AG^AU + PoxtAl/ - AQ. 



(9) 



To further simplify equation @, we substitute for AU us- 
ing the first law of thermodynamics for a system in which 
any work done by gravity is immediately thermalised. In ad- 
dition, even though there is no mass fiow across the outer 
boundary of the system, we still include a mass transfer term 
(fiAN) to account for cooling gas that drops out of the hot 
atmosphere to join the cool phase described previously. At 
this point is it is important to note that the cool phase is as- 
sumed to act as a reservoir that will return mass to the hot 
phase whenever the global heating rate exceeds the cooling 
rate. 

The general form of the first law of thermodynamics 
can be written 



AU ^AQ~ PextAV - A* + ^lAN, 



(10) 



where jj, is the energy per transferred particle, and A7V is 
the change of number of particles in the system, due to mass 
dropping out of the atmosphere, or being added to it. It is 
important to note again that —A*!' is the change of internal 
energy that arises due to the gravitational work being ther- 
malised. By substituting equation (|10|) in to equation (|9]), 
the change in GFE of the hot atmosphere can be written 



AG = -A^ + fi, 



,AN. 



(11) 



Therefore, as the gas cools and loses gravitational potential 
energy (A^ < 0), the GFE will increase unless mass drops 
out of the atmosphere such that |/iatmos AA/"| > |A5'|. 

To make equation (jlip directly applicable to a galaxy 
atmosphere, we express the gravitational and mass loss 
terms as functions of the global heating and cooling rates. 



In the limit that self-gravity of the atmosphere can be ig- 
nored, the change of gravitational potential energy due to 
mass fiow within the galaxy's darkmatter gravitational field, 
during an interval At, can be written as 



A* f« amtotfToM^t, 



(12) 



where riitot is the mass flow rate within the hot atmosphere, 
(TDM is the velocity dispersion of the darkmatter potential, 
and a is a numerical constant the quantifies the gravitational 
potential energy of the gas within f2 as a fraction of a^M- 
The total mass flow rate is determined by the difference 
between the global heating and cooling rates in the atmo- 
sphere, as shown in equations ((2]) and ([3]). The energy flux 
due to mass transfer between the hot atmosphere and the 
embedded cool phase is derived using the following consid- 
erations. The net particle flux into or out of the cool phase 
in a time interval At will be AA*' = /3mtot Af/(m), where 
/3mtot is the mass flux at which hot gas cools and drops out 
of the atmosphere. For convenience, we have expressed this 
mass dropout rate as a fraction of the bulk mass flow rate, 
rhtot, which occurs whenever the global cooling and heating 
rates are not equal, as seen in equation ([2]). A consequence 
of this argument is that a mass fraction 1 — /? of the hot 
atmosphere does not cool signiflcantly, despite being part of 
the bulk flow. Gas which drops out of the flow is assumed to 
join the cool phase. As mentioned above, this acts as a reser- 
voir that can return material back to the hot phase when the 
global heating rate exceeds the cooling rate. For simplicity, 
we assume that the mass returned from the cool reservoir 
makes up the same fraction, /3, of the bulk mass flow rate, 
rhtot. Consequently, the expression AN — l3mtotAt/{m) ap- 
plies whether the gas is being cooled or heated. 

The particles are assumed to be removed from and re- 
turned to the flow at the temperature of the hot phase, 
meaning that /iatmos ~ 'jkhT/l'y — 1). For completeness, a 
net heating rate of j3[H — L) is required to re- heat cool phase 
material up to the temperature of the hot atmosphere at a 
mass flow rate of /Jrhtot. This represents the dominant in- 
put of energy required to return cool phase material back 
to the hot phase flow, i.e. since the hot phase flow is highly 
subsonic, the bulk kinetic energy of the heated material is 
negligible in comparison to its thermal energy. Then, ac- 
counting the work done by gravity and mass transfer, the 
change in GFE can be written 



AG = {H- L)I3 



1 I T\f ^OM 



At. 



(13) 



Accordingly, we see that the system will favour strong heat- 
ing (i.e. H > L) ii (a//?)(7 — 1)(7dm > Cg, i.e. if the gas 
temperature is cool compared to the effective dynamical tem- 
perature of the gravitational potential. In this case, the in- 
crease of internal energy resulting from the work done by 
gravity exceeds the energy that is lost by mass dropping out 
of the atmosphere. In contrast, the system will favour cool- 
ing (i.e. L > H) ii {a/l3){j — l)a^M > c^, i.e. if the gas 
temperature is warmer than the effective dynamical temper- 
ature of the gravitational potential. In this case, increase of 
gas internal energy resulting from the work done by gravity 
is less than the energy that is lost by mass dropping out of 
the atmosphere. 

Assuming that the heating is provided by an AGN at 
the centre of the galaxy then, provided the black hole is 
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massive enough, strong heating (i.e. H > L) will occur at a 
lower accretion rate for radiatively inefficient accretion than 
for radiatively efficient accretion. As described earlier, this 
is because the kinetic outflows couple much more strongly 
to the ambient gas than the photons produced by radia- 
tively efficient accretion. Since the black hole accretion rate 
is likely to increase from a very low value during a cooling 
phase (i.e. L > H), this suggests that the first incidence 
of strong heating will occur at a low Eddington fraction, 
meaning that the AGN power output will be in the form 
of kinetic outfiows. While H > L, the black hole accretion 
rate will be prevented from increasing significantly thereby 
ensuring that the AGN persists in the radiatively inefficient 
state. If this argument is correct, it suggests that internally- 
driven AGN outbursts should be predominantly in the form 
of kinetic outflows. However, in the case that radiatively 
inefficient accretion is insufficient to provide strong heating, 
the accretion rate will continue to increase, causing the AGN 
to enter the radiatively inefficient accretion state. Since the 
power output from the AGN in this state is less effective at 
heating the ambient gas, the accretion rate will continue to 
increase. As a result, the black hole will grow rapidly until 
it is massive enough for radiatively inefficie nt power output 
to pro vide strong heating, as described in IChurazov et al.l 
ll2005ll . 

Importantly, differences between the global heating and 
cooling rates will modify the value of (a//J)(7— l)(crDM/cs)^. 
This is because a, being a function of the gas density 
within Q, Cs and possibly will depend on the historical bal- 
ance between heating and cooling, as well as the under- 
lying gravitational potential. Therefore, we should expect 
(a//3)(7 — 1){(Jy)m/cb)^ to decrease when the gas is heated 
overall, and to increase when the gas is cooled overall. Con- 
sequently, the energetics of the atmosphere favour a heating 
episode while (a//3)(7 — l)crDM > c^. However, when the 
gas has been heated sufficiently that (a//3)(7 — 1)o-^m > , 
the energetics will favour a period during which cooling is 
favoured. This state will persist until {a/l3){'y — 1)o-^m < c^, 
at which point another heating episode will be favoured. 

Equation (|13|l suggests a general solution to the char- 
acteristic differences of AGN activity observed in co ol-core 
(CC) and non cool-core (NCC) galaxy clusters (e.g. iBurnsI 
I1990I : iMittal et al.|[2009h . According to equation (|T3)), the 
lower fraction of AGN outbursts in NCC clusters compared 
to similar CC clusters can be explained by differences in 
the ratio of the dynamical temperature of the gravitational 
potential and the thermal temperature of the atmosphere. 
Therefore, we should expect that NCC clusters spend more 
time with (a//3)(7 — 1)((Tdm/cs)^ < 1, while CC clusters 
spend more time with (a//3)(7— 1)((tdm/cs)^ > 1. These dif- 
ferences could arise due to the formation histories of galaxy 
clusters; for example, NCC clusters may have experienced 
more energetic mergers than CC clusters, thereby disrupting 
the cooling flow and increasing the thermal energy of the gas 
relative to the dynamical te mperature of the gravitational 



potential (e.g. iMcGlvnn fc Fa bian f984; Burns ct al. i997 



iGomez et al.ii200^: iPoole et aLli2008l ). The effect of thermal 



conduction (e.j 



,, Guo et al.l 20081) and differences in metal 
content (e.g. lDubois et aL 20091 ) mav also play a role in this 
behaviour. 

However, while equation (|13p describes the change of 



energetics of gas in the galaxy. In particular, due to its lower 
thermal velocity, the cool phase material can be gravitation- 
ally captured by stars and the supermassive black hole. As 
described in the introduction, the gravitational binding en- 
ergy of gas captured by a compact object is a non-negligible 
fraction of its rest mass energy. Since this is much greater 
than its binding energy to the darkmatter gravitational field, 
some of the cooling gas can effectively gain energy. Then, as 
the gas falls onto the compact object, the binding energy 
will be radiated as heat, thereby reducing the net cooling 
rate of the atmosphere. It is important to note here that the 
total potential energy does not increase, it only appears to 
because the rest mass energy of the gas becomes accessible 
in the presence of strong gravitational fields. 

To account for the influence of the stellar and black 
hole galaxy components, we return to the composite galaxy 
model which describes the 3 main gravitational influences on 
the gas: i) the darkmatter gravitational potential; ii) the stel- 
lar component; iii) the supermassive black hole. The galaxy's 
gas is assumed to move between each of these components 
as the system evolves towards a thermodynamic equilibrium 
state. Therefore, since these 3 galaxy components are closely 
coupled by mass and energy transfer, the overall system will 
tend towards equilibrium as long as the sum of the compo- 
nent free energies decreases. 



2.3 Gibbs free energy of the composite galaxy 
system 

In this extension of the arguments presented above, we make 
the approximation that the cool phase gas will be gravita- 
tionally captured by stars and the supermassive black hole. 
Consequently, each of the 3 gravitational components will 
contain its own quantity of gas, as shown in figure [T] Thus 
there is a GFE for each gravitational component. Further- 
more, due to cooling and heating processes, mass and en- 
ergy will be re-distributed between the components, mean- 
ing that the component GFEs will change as the system 
evolves towards thermodynamic equilibrium. It is important 
to note here that the vast majority of the cool phase gas is 
assumed to be close to infinity with respect to the black 
hole and the stars. As a result, it takes a negligible quan- 
tity of energy to unbind the gas and re-distribute it among 
the different gravitational components. Then, since the gas 
components are closely coupled, the system will tend towards 
equilibrium if the sum of the changes of the component GFEs 
IS less than zero, i.e. J]] AG < 0, even if the GFE of any 
component increases. 

The GFE of the stellar and black hole components each 
take the same form and consist only of mass transfer terms 
which represent the change of energy as mass flows between 
the components. Accordingly, the change of GFE for gas 
that becomes bound to stellar gravitational fields is 



AG* = M*AiV. 



(14) 



where /x« is the energy per particle of gas that is bound 
to the stars, and AA^, is the particle flux into the stellar 
component. The energy gain due to thermally cool gas be- 
ing captured by stellar mass gravitational fields, during an 
interval At, can be written as 



GFE for the hot atmosphere, it does not fully represent the 



li*ANt = ri„r),iht,c At, 



(15) 
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Figure 1. Illustration of one hemisphere of the 3-component galaxy model consisting of the atmosphere, stars, and the black hole. As 
the atmosphere cools, mass will flow inwards, with a fraction dropping out to join the cool phase (indicated by the clouds) which can 
be captured by the stars and the supermassive black hole. The boundary between the atmosphere and the environment is at a pressure 

Poxt. 



where rh* is the rate at which mass is captured by stellar 
mass gravitational fields, and 77* represents the average 
gravitational binding energy of gas to this component as a 
fraction of its rest mass energy, with c being the speed of 
light, and rjc* is the fraction of energy released that could 
couple to the ambient gas. 

The change of GFE for cool gas that becomes gravita- 
tionally bound to the supermassive black hole is 



AGbh = hbh^Nbh, 



(16) 



where ^bh is the energy per particle of gas that is bound 
to the black hole, and AA^bh is the particle flux into the 
back hole component during a time interval At. As above, 
the gain in energy of gas being captured by the black hole 
during the interval At, can be written as 



f^BuANBH = VcVBHrilBHC At, 



(17) 



where ttibh is the rate at which mass is captured by the black 
hole gravitational field, and 7?bhc'^ represents the average 
gravitational binding energy of gas to this component as a 
fraction of its rest mass energy, and rjc is the fraction of 
energy released that could couple to the ambient gas. 

In this model we assume that all of the mass entering 
the cool phase is captured either by the stars or the black 
hole, therefore it follows that Prhtot +m* +tobh = 0. Then, 
if the mass joining the stellar component is a fraction (1 — /) 
of the mass entering the cool phase, we have m, = — (1 — 
/)/3mtot; correspondingly, the mass captured by the black 
hole must be a fraction, /, of the mass entering the cool 
phase, rriBH = -fPrhtot- 

Combining the relations above means that the sum of 
the component GFEs is 



— arhtoto-oM + /3mto 



(7 



+ 



r7c.r;*(l — /)/3mtotC - r?c7;BH./'/3mtotC 



At 



At. 



(18) 



Using the definition of mtot from equations ^ and the 
total change of the GFE is expressed in equations (|19l ) and 
H20p . Specifically, these equations give the change in GFE for 



gas which is being heated at a total rate H, cooled at a rate 
L, accounting for the thermalisation of gravitational work 
and mass transfer to/from the atmosphere, and the gravi- 
tational infiuence stars and the supermassive black hole. 



l-(^ 

Cs 



At, 



where 



(7-1) 



2 2 2 

■gf^DM + (1 - ,f)'7c*7?.C +/j7c??BHC 



(19) 



(20) 



Here, v'^ can be considered to be the energy per unit mass 
gained by the system due to the processes associated with 
gas cooling. Under the assumption that most internally- 
driven accretion events produce kinetic AGN outflows which 
dissipate their energy within the region of interest, the best 
choice for j^c is ---^ 1. Furthermore, the value of 77,2^ in equa- 
tion (|20p is equivalent to the kinetic energy per unit mass 
of stellar winds or supernov a ejecta used in semi-analytic 
galax y evolution models (e.g. ICroton et aLlbOOSlTDave et al.l 
[2OO8). Comparison with such models suggests that (1 — /)??* 
will be much smaller than /?7c7?bh, provided / is greater than 
~ 10"^, meaning that the value of r^c* is generally less im- 
portant than rjc. 

Since the system favours changes that reduce the total 
GFE (J2 AG < 0), heating wiU be favoured {H > L) ii v > 
Cs. That is, heating is favoured when the gain of energy that 
occurs when gas is captured by stars and the black hole -|- 
the thermalisation of gravitational work, exceeds the loss of 
energy that arises when gas cools out of the hot atmosphere. 
The system can also favour cooling (L > H), but only if w < 
Cs. That is, cooling will be favoured if the build up of energy 
in the presence of the stars and the supermassive black hole 
-I- the thermalisation of gravitational work, is less than the 
loss that arises when gas cools out of the atmosphere. It is 
worth noting here that if ?7c ~ 0.05 the system is more hkely, 
for a given value of /, to be in a cooling state since /??c??bhc^ 
will be smaller relative to Cg. Figure [2] shows the regions of 
parameter space which are favoured by systems that can be 
described using equation (|19p . 

For the reasons outlined in the previous subsection, it 
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is important to remember that differences between heating 
and cooling rates will modify the value of (v/cb)^ such that 
we expect [v/csf' to decrease when the gas is heated over- 
all, and to increase when the gas is cooled overall. Therefore, 
the energetics of the system favour a heating episode while 
> Cg. However, when the gas has been heated sufficiently 
that > Cs , the energetics will favour a period during which 
cooling is favoured. Again, this will occur until v'^ < cf, at 
which point another heating episode will be favoured. This 
means that the time-averaged heating rate will closely match 
the time-averaged cooling rate. Consequently, the system 
will evolve quasi-statically towards a thermodynamic equi- 
librium state. Quasi-static evolution has the observable ap- 
pearance of feedback heating regulating the growth of its host 



3 DISCUSSION 



Equation (|19|) is a general expression that describes the cir- 
cumstances under which feedback heating will be favoured. 
However, although the formalism predicts that there should 
be periods of overall heating and cooling, it does not make 
any explicit predictions about the duration or frequency of 
these episodes. In order to solve this problem it would be 
necessary to derive and solve a differential equation for the 
time-evolution of u^, as defined in equation (|20p . Any dif- 
ferential equation of this type wo uld be highl y uncertain. 

Indeed, it was for this reason that iPopd ii2oiir) used a dif- 
ferent approach: to find a Lagrangian function that yielded 
periodic episodes of feedback heating th at were in approx i- 
mate agreement with the observa tions of lBest et al.l l|2005l ): 
iBestI (120071 ) : IShabala et all (120081 ). Notably, the Lagrangian 
that provided the best agreement with observations favoured 
minimum energy output from the system which, in turn, 
favours a scenario in which time-averaged heating balances 
cooling. Below, we demonstrate that the minimum energy 
output Lagrangian actually provides an alternative way of 
expressing the Gibbs thermodynamic potential for this type 
of system. 

Using equation (|19p , we can exp ress the minimum out- 
put energy integral from |Pop3 l|201ll ) in terms of the GFE, 



E,^t{t)^ / [H{t') + L{t')\At' 



(21) 



where G(t) is the rate of change of total GFE with re- 
spect to time. For a system that is evolving steadily and 
quasi-statically, we assume that G = constant < 0, while 
[1 — (u/cs)^] changes sign as a result of i) the work done 
by gravity and ii) feedback, such that over long times (and 
many sign switches) G/\l — {v/csY] At' « 0. In this limit, 
equation ((2T|) reduces to 

\H{t') + L{t')]At' ^ f2L{t')At', (22) 
'0 Jo 

meaning that that time-averaged heating and cooling rates 
are approximately equal: (ff) ~ (L). Thus, the minimum 
energy output from a galaxy is attained when the time- 
averaged heating balances cooling such that the system 
evolves quasi-statically towards equilibrium. This finding 



can explain the seemingly close match between heating and 
cooling rates in elliptical gal axies and galaxy clusters in- 
ferred from observations (e.g. iBirzan et ah Il20q4 iBest et al.l 
l2006l: iNulsen et al.ll2007l : iDunn fc Fabianll2008l ). 

iPopd (|201ll )" showed that a direct consequence of the 
minimum energy output condition is that optimal heat- 
ing occurs in discrete events at periodic i ntervals, in qual- 
itative agreement with observa tions (e.g. iBest et ah I I2OO5I : 
lBestll2007l : IShabala et al.l [20081 ). While this is a useful re- 
sult, it is important to remember that it is an idealised 
heati ng scenario: in reality AGN duty c ycles will be noisy 
(e.g. iKing fc Priiigi^ l20(36l : iPopd l2007l : iPavlovski fc Popd 
I2OO9I : iGaspari et al.ll20l'll )" due to the impact of small-scale 
processes, but displaying (to varying degrees) the features 



ISilk & Reedll998l: King||2009l: 


Pope et al. 201 ll). Bondi ac- 


cretion fe.g. Allen et al.l 2006 


: Cattaneo & Tevssieij l2007l: 


iBooth & Schave 20091; McCarthy et al.l|201C 


), or cold feed- 
) subject to the 


back fe.glPizzolato & Sokeij|2005l: Sokedbooe 


optimality criterion. 



3.1 



A comment on black hole accretion and 
entropy 



In any discussion about a system evolving towards thermo- 
dynamic equilibrium, it is informative to consider entropy 
changes. For a galaxy, it can be shown that the entropy 
increase associated with black hole growth far exceeds any 
entropy losses that occur as a result of gas cooling, or gains 
that result from mass capture by stars. Specifically, the en- 
tropy change when a black hole grows by an energy mc? 
is ( e.g. lBekensteinlll973l : iBardeen et al.lll973l : lHawking||l976D 



A5b 



mc^At 
Tbh ' 



(23) 



where Tbh is the black hole temperature given by (e.g. 
lHawkinel[l97^ 



Tbh 



hc-^ 



(24) 



87rGNA/BHfeb 

with h being the reduced Planck's constant, Gn is Newton's 
gravitational constant, and Mbh is the black hole mass. It 
can be seen from equation (|24|l that the black hole tem- 
perature falls as its mass increases, thereby increasing the 
entropy gain whenever the black hole accretes material. 

In the limit that the black hole dominates the total 
heating rate, H, the power output and black hole accretion 
rate are related hy H = 7;c'7BH?7^BHC^. This means an energy 
(1 — ??BH)rhBHC^ = (1 - riBii)H/{r]cr]Bii) is captured by the 
black hole. Therefore, equation (|23p becomes 

(1 -77bh) HAt 



ASb 



(25) 



VcVbh Tbh 

Then, since the temperature of the black hole is much lower 
than that of the hot atmosphere, accretion onto the black 
hole causes the entropy of the galaxy to increase while the 
GFE decreases. 



4 SUMMARY 

The aim of this article has been to better understand the 
origin of feedback heating processes that are fundamental 
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(H-L)/L 



Disfavoured 



Favou red 
(cooling state) 



Favoured 
(heating state) 



Figure 2. Parameter space for which the Gibbs Free Energy decreases. Systems experiencing feedback heating wiU occupy the shaded 
region in the top right of the figure, where H > L, and (v/cg)^ > 1. In contrast, systems that are cooling will occupy the shaded region 
in bottom left of the figure, where L > H, and {v/cs)^ < 1. A system will switch between these favoured states, depending on the ratio 
(v/cb)^ which will vary over time as the gas is heated and cooled. The y-axis ranges from -1, for H = 0, to an arbitrary upper limit of 
+4, which implies H = 5L. Similarly, the a;-axis ranges from 0, for v = 0, to an arbitrary upper limit o{ {v/csf = 3. The area of a region 
is not an indication of the probability of the system occupying that state. 



in governing the evolution of gas in massive galaxies and 
galaxy clusters. To achieve this, the Gibbs thermodynamic 
potential was calculated for a generalised galaxy model that 
describes the 3 main gravitational influences on the galaxy's 
gas: i) the dark matter gravitational potential; ii) the stars; 
iii) the supermassive black hole, as shown in flgure[T] Differ- 
ences between global heating and cooling rates result in the 
movement of gas within the galaxy, re-distributing mass and 
energy between the 3 gravitational components. Using this 
approach, it was shown that periods of heating and cooling 
are required to ensure that the re-distribution reduces the 
free energy of the gas overall, such that the system evolves 
always towards thermodynamic equilibrium. The key details 
and flndings are described below. 

(i) As the hot atmosphere of a massive galaxy cools by 
radiating X-rays, it will lose gravitational potential energy 
while simultaneously being compressed by the overlaying 
gas. In addition, the cool gas will be captured by the stars 
and the supermassive black hole at the centre of the galaxy. 
All of these processes will increase the energy of the gas, 
while cooling gas that drops out of the atmosphere will re- 
duce the thermal energy. To ensure that the GFE always 
decreases, strong feedback heating will be favoured when- 
ever the rate at which the gas gains energy (through cool- 
ing) exceeds the rate of energy loss that occurs when gas 
drops out of the atmosphere, i.e. when (v/cb)^ > 1. Under 
these circumstances, the galaxy will be in a heating state, 
as illustrated by the top right panel of figure [S] 

(ii) When the gas is heated overall, we expect v/cb to 
decrease. Therefore, at the point that (v/cb)^ < 1, the 



system will enter a period during which overall cooling is 
favoured. This is because the rate at which the gas gains 
energy (through cooling) is less than the energy loss that 
occurs when gas drops out of the atmosphere. Under these 
conditions, the galaxy will be in a cooling state, as shown in 
the bottom left panel of figure[2] Furthermore, while the gas 
cools overall, we expect v/cb to increase. Therefore, at the 
point that (v/cb)^ > 1, the system will enter another period 
during which overall heating is favoured, and so on. 

(iii) As a direct consequence of the above, the energet- 
ics of the system favour a close balance between aver- 
age heating and cooling rates, explaining correlations be- 
tween inferred heating and c ooling rates i n elliptical galax- 
ies a r id galaxy clusters (e.g. iBirzan et al.ll2004l : iBest et aP 
l2006l : lNulsen et al.ll2007l : lDunn fc FabianlbOOsf KMore signif- 
icantly, this means that the galaxy evolves quasi-statically 
towards a state of thermodynamic equilibrium. It is for this 
reason that AGN feedback regulates the growth of massive 
galaxies. 

(iv) It has also been shown here that minimising the GFE 
of gas in a galaxy is consistent with minimising the sum 
of the heating and cooling outputs from the system, thus 
providing a broader theoretical mo tivation for the optimal 
AGN heating scenario proposed by IPopd (j201ll l. 

(v) Cycles of AGN activity can be explained by galax- 
ies flipping between the heating and cooling states shown in 
figure [2l As described above, this flipping is driven by the 
variation of the parameter which is influenced by global 
heat ing and cooling rates. By comparison with observ ations 
(e.g. iBest et ahllioosl : lBestll2007l : IShabala et al.1120081 ). mas- 
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sive galaxies with larger radio AGN duty cycles must spend 
more time in the heating state than lower mass galaxies. 

(vi) The same effect can explain the characteristic dif- 
ferences of AGN feedback in cool -core (CC) and non cool- 
core (NCC) galaxy clusters ("e.g. lBurnj|l990l : iMittal et all 
|2009( ). According to this model, a central AGN located in 
a NCC clusters must spend more time in the cooling state 
than a central AGN in a CC cluster. This could occur be- 
cause of different formation histories meaning that NCC 
clusters have systematically higher gas temperatures rela- 
tive to the dynamical temperature of the cluster gravita- 
tional pot ential, or becaus e of the effect of thermal conduc- 
tion fe.g. IGuo et al.ll2008l ). or differences in metal content 
(e.g. [Dubois et al.ll201ll ). 
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